Non-equilibrium Josephson— like effects in mesoscopic S-N-S junctions 
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Wide mesoscopic superconducting-normal-metal-superconducting (S-N-S) junctions exhibit Andreev 
bound states which carry substantial supercurrents, even at temperatures for which the equilibrium 
Josephson effect is exponentially small — the currents carried by different states can cancel each other. 
This cancellation is incomplete whenever the junctions are driven out of equilibrium, e.g., by a dc 
voltage. This leads to phenomena similar to the usual dc and ac Josephson effects, but dominated by 
the second harmonic of the Josephson frequency, which may explain some striking recent experiments. 
A simple description of these, in the spirit of the Resistively-Shunted-Junction model, is suggested. 

PACS numbers: 74.50.+r, 74.40. +k, 74.80.Fp, 73.23.Ps 
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Mesoscopic superconducting-normal- metal-supercon- 
ducting (S-N-S) junctions exhibit a spectrum of low- 
lying electronic "Andreev bound states" , which depends 
strongly on the phase-difference <p between the order pa- 
rameters in S (see Fig. [I]). This ^-dependence persists 
when the temperature T is raised, and the "normal- 
metal coherence length", = ^hD/ksT, becomes 
much smaller than the distance L between the two S elec- 
trodes (D is the diffusion constant in N). In this high- 
temperature regime the equilibrium Josephson coupling 
is exponentially weak. However, as the spectrum depends 
on <j), and hence on time t (the voltage V is proportional 
to dtfi/dt ), non-equilibrium (NEQ) effects occur, includ- 
ing supercurrents which are the topic of this Letter. 

Such junctions have attracted much attention. The </>- 
dependence of their conductance was studied extensively 
in recent years jj| , and it was demonstrated that a naive 
Ginzburg-Landau description of proximity effects was in- 
sufficient. In these "Andreev interferometers" N is con- 
nected to external normal-metal electrodes, which ther- 
malize the occupations 0. In contrast, for an "isolated" 
N the electronic excitations are "bound" : an electron 
near the Fermi level in N cannot enter S because of the 
superconducting gap, except by "Andreev reflection" - 
it pairs with another electron from the Fermi sea, leaving 
behind a hole. Conversely, a hole may break up a Cooper 
pair and produce an electron. These processes coherently 
mix the electron and hole states in N, in a ^-dependent 
manner. As energy-relaxation is slow in the mesoscopic 
regime (t p > L? / D, with t v the single-particle dephas- 
ing time), NEQ situations naturally develop. 

NEQ phenomena have been studied for both clean 
|,| and dirty |,§ S-N-S junctions §. The progress 
made here is in (a) developing a simple but versatile 
model of NEQ phenomena, in the spirit of the resistively- 
shunted-junction (RSJ) model (at the price of a restric- 
tion to small voltages); and (b) using the unusual har- 
monic content of the NEQ supercurrents to understand 
a recent observation , which cannot be explained M in 




FIG.g: Energies of some of the Andreev bound states E n 
(equally spaced in n), as a function of <j), for a diffusive 
S-N-S junction (see inset) with 5 <C Ec <C A (and t v — > 
oo). Thick line: a representative used below, E rcp ((j>). 

terms of previously proposed effects fic|j . The simplest 
description of NEQ, with a constant energy-relaxation 
rate 1/te (with te > T v ), is used to discuss wide (or 
long) dirty junctions, in the high-temperature regime, 
E c = fiD/L 2 « T « A. Here E c is the Thouless en- 
ergy (the correlation energy of the spectrum), ks = 1, 
and A is the gap in S. The single-particle level spacing 
in N, S, is taken small (metallic limit): 5 -C Ji/te <S Ec- 
The density of states u(e, <j>) can be calculated from 
the (disorder-averaged) Usadel equations |Tl]Jl^ ]. It is 
convenient to define the energy E n (cj>) of the nth level: 
n = Jq™ f(e, 4>) de, see Fig. [J (e = is the Fermi level). 
We have assumed for simplicity that time-reversal sym- 
metry is preserved, that the pair potential A is constant 
in S and vanishes in N, and that the N-S boundaries 
are "perfect" , i.e. that no normal reflection occurs and 
the amplitude for Andreev reflections is of unit magni- 
tude. A mini-gap is proximity-induced in N. Its size is 
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E g ~ Z.lEc for 



0, and it closes when 



7r and 



reopens periodically [[L2| (fluctuations in E g and -En are 
only ~ (5). The spectrum is very sensitive to symmetry 
breaking — for example, a moderate amount of spin-flip 
scattering, with h/r s f ~ -Ec> can close the mini-gap. 
However, for most S-N-S geometries it contains a band 
of width ~ Ec of levels with a significant (^-dependence. 
For our system, the curve E Iep ~ 3AEc \/ 1 + 0.7 cos(</>) 
represents the low-lying part of the spectrum well, and 
the effective number of levels in this band is TV ~ 10Ec / S, 
including a factor of 2 for spin [a simple shift in en- 
ergy does not affect the contribution of a level E n (<j>)]. 
As we will use only the fact that E Yep is </>-periodic, 
the results will be qualitatively applicable also to the 
experimentally-relevant clean case, with a mean-free- 
path of a few times L (though perhaps not to the theo- 
retical clean limit ||Q, which has a separable spectrum). 

The states E n carry currents, I — (2e/h)(dE n /d(j)) per 
occupied state, as seen by equating the spent and stored 
energies, IV dt = Yl dE n (neglecting any changes in the 
interaction energy), and using the Josephson relationship 
d<f> = (2e/K)Vdt. For each n there is also a state at — E n , 
with the opposite current. The total supercurrent is thus 



(1) 



where / is the occupation probability of E n , and 1 — / is 
that of —E n (Ref. the spin index is included in n). 

The fact that NEQ occupations often enhance super- 
currents was demonstrated in the seventies ||r), and fol- 
lows from an elegant argument. In thermal equilibrium, 
/ = / cq = 1/ (l+exp(e/T)), and Ig takes the form 



I eq (cf>) = / dej{e,<t>) tanh(e/2T) , 



(2) 



where j(e) is the "Josephson current density" , j(E n ,4>) oc 
v(E n , <p)(dE n /d<f)). This j(e) is the imaginary part of a 
"Green's function" which is analytic in the upper half of 
the complex e plane — a positive imaginary part of e cor- 
responds to a dephasing rate (h/r^,) which would smooth 
out any singularities. Using contour integration (the in- 
tegrand is even), one finds the well-known Matsubara 
sum: 7 oq = 2iriTJ2T j(i u, n)t where ui n = (2n— l)irT 
are the poles of the tanh factor. At high temperatures, 
even the smallest Matsubara frequency has a decay time 
h/u>i < L 2 /D, yielding an exponentially small j(iu} n ). 
Thus, for physical quantities of this form, thermal aver- 
aging is mathematically equivalent to dephasing. 

In NEQ situations contour integration cannot be used, 
and the currents are not exponentially small. For Eq. (|^) 
to give small results, j(e) must oscillate: in our diffusive 
system, j(iuj) decays exponentially with y/cJ/Ec on the 
imaginary axis, and correspondingly j(e) oscillates and 
decays rapidly with y/e/Ec- The E n (<fi) curves in Fig. |l| 
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FIG. ^: (a) Ig vs. t, and (b) amplitudes Ik of its harmonic 
decomposition at t 3> te (with fc > 0; phases not shown), 
for five dc voltages labeled by v = (2e /K)Vte ■ 

thus change their character repeatedly at higher e, occa- 
sionally having shallow maxima at <fi — n, rather than 
minima. As can be seen from Eq. (^) below, the devi- 
ations of / from equilibrium change sign in rhyme with 
these oscillations, and thus the integrand of the NEQ part 
of the supercurrent, 2 J Q dej (/eq — /), does not change 
sign and cannot be affected by any cancellations. 
The expression for Is may thus be approximated by 



Is 



(3) 



The errors incurred here ]l4|| are probably smaller than 
those of the relaxation-time approximation, 



dt te 



/eq), 



(4) 



which we shall also employ. For example, the latter ig- 
nores the effects of the mini-gap on the electron-electron 
and electron-phonon relaxation processes. 

When a dc voltage is applied to the junction, dtfi/dt = 
2eV/h — const., this model gives ac supercurrents, see 
Fig. H (the units 7 ncq and tj are defined below). Here 
T > E c was used to take / cq = and / cq = ^—E rep /4T. 
The second harmonic dominates because Is of Eq. (0) is 
a product of two oscillatory factors, dE rcp /d(j) and (/ — 
/ cq ). For junctions with t v ~ L 2 / D one expects a much 
"softer" spectrum, with i^-cp oc cos(</>), and Is would 
exhibit true frequency doubling, evolving from a sin 2 (0) 
behavior to — sin(2</>) as V is increased. 

In the limit of small voltages, we reproduce the results 
of Refs. §0: /-/ cq ~ -T E (2eV/n)(df cq /d ( f > ) 1 and 
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FIG. 0: (a) The "free-energy" function F(<£,/), for / = 
0, is a "winding valley" with minimum at / = E(<f) and 
a parabolic cross-section, (b) The probability density in 
/, for t e — > oo, is p(f) oc / d(j> exp(—F/T^); here = 
0.25S ncq . (c) The effective potential, F(f) = -T^Xogp. 
When T <CT^, its minima "trap" the system. 



This adds a phase-dependent term to the ohmic conduc- 
tance of the junction. For a large dc voltage, we may 
approximate / by a constant — the phase average of / e q 
— which yields a purely oscillatory Is- As can be seen 
from the "Debye mechanism" of Ref. H , which monitors 
the steady-state transfer of energy into the heat bath 
(oc 1/te) rather than that into the electrons, the total 
power dissipated approaches a constant P m ax at large V, 
and so f I dt decreases as P m ax/V. 

So far, we have ignored the normal current compo- 
nent, In ~ GnV, where Gn is the conductance of N in 
the absence of proximity effects. Clearly, in the metallic 
limit the occupations can not evolve adiabatically (as in 
Ref. ||). Instead / diffuses in energy with a coef- 
ficient De{c,4>) ~ GnV 2 S. To see this, note that in a 
time t a particle can cross the junction ~ tEc/h times, 
shifting in energy by ±eV at each Andreev reflection, so 
that D E ~ (eV) 2 D/L 2 [recall that G N ~ (e 2 /h)E c /S\. 
This additional relaxation term can be neglected, com- 
pared to the terms kept in Eq. (0), if teDe -C Eq or 
tjE*E <C Eq, depending on the context [tj = h/2eV is 
the Josephson period). This limits our model to small 
voltages, eV <C \J Ec^Ite, or eV <C Ec, respectively. It 
is difficult to avoid the latter restriction, because at still 
higher voltages the temporal and spatial dependencies of 
the relevant Green's functions become intertwined. 

In many situations I(t) rather than V(t) is given. 
As in the well-known RSJ model (jl6), we write / = 
Is + In + If, where Ip is a fluctuating noise compo- 
nent (S-N-S junctions are generally overdamped, with no 
displacement current). We rescale E rop (<fi) and /, giving 
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where E and / are defined by E Icp (4>) — AE + B 
and / = | — [Af + B)/4T, with the requirement 



FIG. H The I — V curves and differential resistances (in- 
sets) for a dc current bias in the NEQ model, for four 
values of the relaxation time te', compared to the RSJ 
model with I c = I ncq , and to I = GnV (dotted lines). 

— 1 < E < 1 fixing A and B (in our model A ~ 
1.3E C ). Here I ncq = (2e/h)NA 2 /4T 7 and the energy 
scale is -Encq = (fr/2e)/ noq — 4.1Eq/ST, only a factor 
of Ec/3-5T smaller than the T — Josephson coupling 
energy. The correlators (I F (t)I F (0)) = 2T 4> G 8{t) and 
(Jj?(t) Jf(0)) = (t e T / E neq ) 5(t) give the Gaussian fluc- 
tuations of Ip and Jp. For thermal noise, — T; high- 
frequency external noise can be described by > T . 

Eqs. (0) and (Q) describe overdamped motion in a free- 
energy |Lp landscape F(<j>, f) = E ncq (E-f) 2 - (h/2e)Icf), 
see Fig.T3|(a). This generalizes the tilted washboard po- 
tential of the RSJ model, F(<j>) = -(h/2e)(I c cos cA + 1(f)) 
(where / need not be followed). For fixed <f>, the fluctu- 
ations induced by Jp reproduce the y/N noise of ther- 
mal occupations (and result in increased current noise, 
I s oc dF/dcj), Ref. @). Note that for t e -> (fast equi- 
libration), / may be integrated out, which must repro- 
duce the RSJ model. Here I c = 0, but the ^-dependence 
of the "bottom of the valley" re-emerges at lower T. 

In the opposite limit, te — > oo, the "fast variable" <fi 
may be removed (see Fig. ||), giving a two-peaked prob- 
ability density in / (one broad peak for large T^). The 
corresponding effective potential, F(f), has valleys with 
a depth roughly ~ T^, which trap the system if T << 
(Ref. This drastically reduces the rate of phase- 

slip: the <fi variable remains near the extrema of E, and 
changing it with / fixed (because te is large) entails over- 
coming a barrier of order E ncq . This surprising fact - 
that noise enhances conductance near I — — is demon- 
strated numerically below. 

The dc I — V curves of this model, obtained with 
I(t) — const, and no noise (Ip = Jp = 0), are dis- 
played in Fig. [| using units of J neq , V neq = G N I neq and 
Tneq = Ti/2eV neq for current, voltage and time. The re- 
sults for large Tp lie remarkably close to the I — V curves 
of the RSJ model (^0|, except for the finite zero-bias 
conductance, equal to 1 + r, or Gat + TE{2e/h)I ncq . Re- 
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FIG. |: (a) The effect of external noise (indicated by 
T<f,) on the dc / — V curves with r = 25 (the waviness is 
due to the finite numerical integration time). Inset: on 
a much expanded scale the = 0.25 curve is seen to 
cross the one for — 0, and appears to have a "critical 
current" of ~ 0.1/ noq . Thermal noise counteracts this 
(T > 0; dashed lines), (b) An ac current bias (with 
amplitudes 7 ac , frequency ui — 2eV ncq /h, and no noise) 
produces Shapiro steps. The step at V = 0.5 is due to 
supercurrents at twice the Josephson frequency. 



suits with external noise (and a large te) are shown in 
Fig. ||(a). Such numerical calculations indicate that for 
finite te the scaled conductance grows at most as r 2 , and 
not exponentially as in the te — > oo case. 

An ac drive, I ac cos(wt), produces Shapiro steps in the 
/ — V curves — phase-locking of the ac Josephson oscil- 
lations to the external frequency u>, see Fig. |(b). As m 
the RS J model [|l6| , closely-related results obtain for the 
simpler, voltage-biased case For weak ac driving, 

the sub-harmonic step at hu> dominates. As observed in 
the experiments of Ref. ||, this step is the largest for 
T ~ Ec , but does not appear to be parametrically larger 
than the others (a mechanism with true frequency dou- 
bling would behave differently). The experimental steps 
were much smaller in size, which could indicate fast equi- 
libration — the steps disappear quadratically for small 
te, because to leading order Eq. ^ gives I — G((j>)V, 
which relates the time integrals of the current and volt- 
age by a constant ratio, J G(<f>) d(f)/2ir, and precludes any 
steps || . A detailed modeling of the experiment requires 
additional measurements, which are in progress. 

In summary, we have discussed NEQ effects very sim- 
ilar to both the dc (Fig. ||) and the ac (Fig. ||) Joseph- 
son effects, which decay with temperature only as / neq oc 
Ec/T. The NEQ occupations can be produced by mi- 
crowave irradiation, a noisy external circuit, or simply an 
applied voltage. The new effects have the signature of be- 
ing dominated by the second harmonic of the Josephson 



frequency, and may have already been observed. 
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